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Abstract: Ecological interactions of species and thus their spatial pat¬ 
terns may differ between homogeneous and heterogeneous forests. To 
account for this, techniques of point pattern analysis were implemented 
on mapped locations of tree individuals from two 1-ha tropicalforest 
plots in Vietnam. We analyzed the effect of environmental heterogeneity 
on tree distributions; spatial distribution patterns of dominant species; 
inter-specific associations; and conspecific associations between life 
stages. Our analyses showed that: environmental conditions were homo¬ 
geneous at plot 1 but heterogeneous at plot 2; in both plots, all six domi¬ 
nant species were aggregated at various scales up to 30 m, and tree spe¬ 
cies were aggregated at larger scales in the homogeneous site than in the 
heterogeneous site; attraction between pairs of species was remarkably 
higher at the homogeneous site while negative associations were more 
frequent in the heterogeneous site; some species, H. kurzii, T. ilicifolia 
(homogeneous plot) and D. sylvatica, S. wightianum (heterogeneous plot), 
showed a lack of early life-stage individuals near conspecific adults. 
Moreover, additional clustering of young individuals was independent 
from conspecific adults, except D. sylvatica in both sites. These findings 
are consistent with the Janzen-Connell hypothesis. Overall, habitat het¬ 
erogeneity influences spatial patterns and inter-specific associations of 
the tree species and evidences of self-thinning are shown in most species. 

Keywords: environmental heterogeneity, spatial point pattern analysis, 
tropical evergreen forest, northern Vietnam. 
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Introduction 

A central question in forest ecology is how to understand the 
processes and mechanisms that control species coexistence and 
community structure, especially at various spatial scales. These 
relevant issues have been addressed in numerous studies. Barot 
(2004) highlighted the impact of both exogenous and endoge¬ 
nous factors on the spatial and temporal distributions of tree 
species. Other studies investigated intra- and inter-specific inter¬ 
actions (Callaway and Walker 1997; Bruno et al. 2003), dispersal 
limitation (Hubbell 1979), negative density dependence (Wright 
2002), or habitat preference (Condit et al. 2000). Tilman (2004) 
emphasized that spatial patterns of trees may be explained by 
environmental niche effects and trade-offs among species due to 
dispersal and competitive ability. 

Environmental heterogeneity (such as different soil types, rock 
outcrops, or streams) makes spatial pattern analysis more com¬ 
plicated because it confounds biotic and abiotic effects (Li and 
Reynolds 1995; Wiens 2000). Dispersal limitation is emphasized 
as a potential mechanism for separating species in space and 
reducing competitive exclusion (Seidler and Plotkin 2006). Be¬ 
sides that, a patchy distribution of trees can also be caused by 
habitat preference where demographic processes and limiting 
resources may simultaneously influence spatial patterns (Wagner 
and Fortin 2005; Getzin et al. 2008). Spatial aggregation of a 
species can be induced by limited seed dispersal or patchy habi¬ 
tat conditions or by both factors (Webb and Peart 2000). 

In addition, negative density dependence or self-thinning is a 
prominent mechanism for regulating population dynamics and 
facilitating species coexistence (Wright 2002). This mechanism 
has been considered by a negative density of conspecific distance 
relations in forest dynamics such as recruitment, growth or sur¬ 
vival (Condit et al. 1992; Peters 2003; Uriarte et al. 2004a). Den¬ 
sity comparisons of life history stages ranging from seedlings to 
adults provide evidence of self-thinning. However, self-thinning 
cannot always be observed in species with distinct habitat pref¬ 
erences (McDonald et al. 2003) or in tree fall gaps (Augspurger 
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1984). 

Direct tree-tree interactions are often observed at small scales 
while at large scales, they can be obscured by environmental 
factors (Wiegand et al. 2007). The challenge is to separately 
consider first-order effects (i.e., habitat preference of tree species) 
and second-order effects (i.e., competition or facilitation within 
or between species). Disentangling these effects makes it possi¬ 
ble to find evidence for different hypotheses about species coex¬ 
istence and structuring of plant communities at different spatial 
scales (Callaway and Walker 1997; Wiegand et al. 2007). Inter¬ 
estingly, Getzin et al.(2008) found that plant ecology, in terns of 
plant population dynamics and pattern fonnation, may differ 
between homogeneous and heterogeneous sites, beyond the 
purely statistical effects of heterogeneity. Given that most studies 
focus on homogeneous sites, there is an urgent need to compare 
spatial patterns in homogeneous and heterogeneous forest sites. 

In this study, we analyzed spatial patterns and associations of 
six dominant tree species in the tropical evergreen forest in Babe 
National Park, in northern Vietnam. We investigated four aspects 
of tree distribution: environmental heterogeneity, spatial distri¬ 
bution patterns of species, inter-specific associations, and con- 
specific associations between life stages. We used methods of 
point-pattern analysis and compared patterns to different null 


hypotheses. Significant deviation from null hypotheses was 
evaluated by confidence envelopes , which were constructed 
with Monte Carlo simulations (Stoyan and Stoyan 1994; Diggle 
2003). This approach has shown theoretical and practical advan¬ 
tages in previous studies on spatial-pattern analysis of tree dis¬ 
tribution (Wiegand et al. 2007; Getzin et al. 2008; Martinez et al. 
2010; Wang et al. 2010; Lan et al. 2012). We expect this study to 
contribute to the understanding of spatial patterns and environ¬ 
mental effects on species distributions in the tropical forests. 

Materials and methods 

Study area 

The study was conducted at Babe National Park (NP), in north¬ 
ern Vietnam. It is located approximately 250 km north of Hanoi 
and 75 km north-west of Backan City. This area receives an an¬ 
nual rainfall of 1378 mm and a mean temperature of 22°C and 
mean humidity of 83.3%. Babe natural lake has a surface-water 
area of about 500 ha and is surrounded by straight cliffs of Karst 
Mountains. Two plots were chosen in the core zone of the park 
where human impact is minimal (Fig. 1). 


Conservation zones in Babe national park 


105°34'0"E 105°36'0"E 105"38'0"E 105"40'0"E 105’42'0"E 105"44'0"E 



Fig. 1: Map of Babe NP and location of study plots (Source: the NP) 


The study site is classified as tropical evergreen lowland forest, 
which is one of several rainforest sub-types in the area (Rcfee 
2011). The soil has a thick layer of yellow-brown ferralsol with 
clay to sandy clay particle-size classes. 
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Data collection 

We selected two study sites that differed in geomorphological 
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characteristics such as slope position (hillfoot and hillside); rock 
outcrops (low and high abundance); and slope (Table 1). In each 
site, a 1-ha plot was established and subdivided into 100 10 m X 
10 m subplots. The diameter at breast height (dbh; at 1.3 m 
above the ground); tree coordinates ( x,y ); and species name were 


recorded for all woody plants with dbh >2.5 cm. Stem-mapping 
of individuals was done using a laser distance measurement 
(Leica Disto D5) and compass. All individual trees were divided 
into three size classes: “juvenile” (dbh <5 cm), “sub-adult” (5 cm 
< dbh <10 cm) and “adult” (dbh >10 cm). 


Table 1: Environmental characteristics of the two 1-ha forest plots in Babe NP. 



Area (ha) 

Coordinates 

Elevation (m) 

Slope (°) 

Aspect 

Rock outcrops 

Location 

Plot 1 

1 (lOOxlOOm) 

22°24'567"N, 105°37784" E 

323 

30 

East 

Few 

Hill-foot 

Plot 2 

1 (100 X 100 m) 

22°25'053 ,r N, 105°37'744 "E 

331 

20 

East 

Abundant 

Hill-side 


For most of the analyses, we focused on the five most domi¬ 
nant tree species in each plot. Four of these five species occurred 
at both plots: Diospyros sylvatica (Ebenaceae); Burretiodendron 
hsienmu (Tiliaceae); Hydnocarpus kurrzii (Flacourtiaceae); 
Syzygium wightianum (Myrtaceae).Two species occurred on one 
plot only: Taxotroplris ilicifolia (Moraceae) and Streblus 
macrophyllus (Moraceae). H. kiirrzi, T. ilicifolia and S. macro- 
phyllus are shadetolerant species, which develop well on thick 
and humid soils, while B. hsienmu and S. wightianum are 
shadeintolerant species (FIPI 1996). 

Statistical analysis 

The K-function and pair-correlation function 
The if-function is defined as the expected number of points oc¬ 
curring within distance r of a randomly chosen point and nor¬ 
malized by the intensity X of the point pattern (Ripley 1976; 
Diggle 2003; Wiegand and Moloney 2004; Illian et al. 2008). 
K(r) is the cumulative distribution function up to a given radius r. 
We used the square root transfonnation of the W-function: 

L(r) = (K(r) / 7r) V2 — r (1) 

The univariate pair-correlation function g(r) is related to the 
widely used .^-function. It is a non-cumulative distribution func¬ 
tion and presents the expected number of points in a ring at dis¬ 
tance r from an arbitrary point divided by the intensity of the 
pattern (Stoyan and Stoyan 1994; Wiegand and Moloney 2004). 
It is calculated as: 

g(r) = K /(2m-) for r > 0 (2) 

The parameter g(r) indicates whether a pattern is (1) random 
(complete spatial randomness CSR; g(r) = 1), (2) clumped (g(r) 
> 1) or (3) regular (g(r) < 1) at a given radius r. It can be ex¬ 
tended to describe the relationship between two point patterns. 
The bivariate pair-correlation function g, 2 (r) represents the ex¬ 
pected density of points of pattern 2 at distance r from an arbi¬ 
trary point of pattern 1, divided by the intensity of pattern 2 
(Wiegand and Moloney 2004). Therefore, g\fr) = 1 indicates 
independence (no interaction) between two types of points, and 
gn( r ) >1 indicates that there are, on average, more points of type 
2 at distance r from points of type 1 than expected under inde¬ 


pendence, hence showing attraction between these two types of 
points. Similarly, gn( r ) <1 indicates repulsion between the two 
types of points. 

Null models 

We used a two-step approach as suggested by Dixon (2002) and 
Wiegand and Moloney (2004): (1) compare the estimated L- or 
g- function (equation 1, 2) with the expected values from an ap¬ 
propriate null model, which is potentially useful for the underly¬ 
ing point processes; and (2) hypothesistesting based on approxi¬ 
mate confidence envelopes constructed by Monte Carlo simula¬ 
tions of the stochastic process. 

Complete spatial randomness (CSR): CSR is a commonly 
used null model for univariate analysis. It assumes no interaction 
between the points of a pattern, for example locations of indi¬ 
vidual trees of a given species or size class (Wiegand and 
Moloney 2004). The intensity of trees in the pattern is constant 
over the study plot, meaning that trees can occur at any position 
and are independent from other trees. 

Heterogeneous Poisson process (HP): HP is used if the dis¬ 
tribution of trees is of varying intensity within the pattern. This 
may be caused by obstacles like rock outcrops, streams or dif¬ 
ferences in soil. In such cases, the null model of CSR is not suit¬ 
able for exploration of second-order characteristics (e.g., 
tree-tree interactions). The heterogeneous Poisson process is an 
alternative to account for the first-order moment (i.e., large-scale 
variation in habitat quality). The constant intensity assumed for 
CSR is replaced by an intensity function that varies with tree 
location but preserves the independence of tree positions from 
each other. For simulations with the grid-based Programita soft¬ 
ware, a circular moving window with radius R is placed at a 
provisional point of a random cell (x,y), then a probability den¬ 
sity is computed and retained for this point. The procedure is 
repeated for all N points of the pattern (Wiegand and Moloney 
2004). 

Random labeling : If we want to investigate whether or not 
pattern 1 and pattern 2 are created by the same stochastic process, 
locations of points in the joined pattern (1+2) are fixed and then 
the locations of point pattern 2 are randomized among fixed lo¬ 
cations (Wiegand and Moloney 2004). The random-labeling null 
model tests if a second stochastic process randomly affects the 
individuals of a population or if there is a spatial pattern in this 
second process (Goreaud and Pelissier 2003). Two test statistics, 
gn-gn andg 2 r.? 22 > were examined. While gn-gn shows different 
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intensities of pattern 2 points around pattern 1 points compared 
to pattern 1 points around themselves in circles with radii r, 
g 2 i~g 22 evaluates whether type 2 points tend to be surrounded by 
points of type 1 (g 2 rg 22 > 0) or more by other points of type 2 
(g2l-g22< 0). 

Analysis 1: Environmental heterogeneity 

An initial and fundamental step in spatial pattern analysis of 
plant communities is to identify if environmental conditions on 
the site are homogeneous or heterogeneous. Tree densities may 
vary between locations due to different environmental conditions 
and tree-tree interactions. These two effects are difficult to dis¬ 
tinguish. 

Our study sites were chosen to account for and compare envi¬ 
ronmental effects assuming homogeneity of environmental con¬ 
ditions at plot 1 and heterogeneity at plot 2. Getzin et al. (2008) 
argue that adult trees of all species have exploited all available 
sites. Thus, the spatial pattern of all adults should clearly reveal 
the effect of environmental heterogeneity on tree distribution. 

On the heterogeneous site, we expected the intensity of trees 
to vary at large scales. Therefore, the g-function should not ap¬ 
proach the expected value 1 for scales which are beyond direct 
interactions of trees. The Z-function should strongly increase and 
not converge to 0 for large scales because it is an accumulative 
function of r (Wiegand and Moloney 2004). Comparing results 
from the g- and /.-functions is considered an advantage over 
visual interpretation of habitat heterogeneity. 

Analysis 2: Species distributions 

We used the univariate pair-correlation function to investigate 
spatial patterns of the six dominant species in both plots. Null 
models were CSR for plot 1 and HP for plot 2. A radius R = 30 
m for the circular moving window was chosen to factor out the 
large-scale structure of the pattern in plot 2. Under the HP null 
model, the locations of the focal tree were displaced randomly 
within the neighborhood R. This process eliminated neighbor¬ 
hood patterns at small scales r < R while preserving larger-scale 
patterns (by randomizing the pattern using its intensity). 

We used a displacement band R = 30 m that is similar to typi¬ 
cal scales at which tree-tree interactions have been found 
(Hubbell et al. 2001; Peters 2003; Uriarte et al. 2004b). We 
combined the non-parametric Epanechnikov kernel for the inten¬ 
sity function with a spatial resolution of 1 m. We expected inter¬ 
actions of tree species at small scales and environmental hetero¬ 
geneity at larger scales. 

Analysis 3: Species associations 

To identify inter-specific association of tree species (i.e., repul¬ 
sion or attraction), we perfonned the bivariate analyses with null 
models appropriate for environmental conditions. For the homo¬ 
geneous site, pattern 1 was fixed and pattern 2 was randomized 
using CSR null model (independence null model). For the het¬ 
erogeneous site, we kept the location of trees in the first pattern 
fixed and randomized the location of trees in the second pattern 
using a HP null model. As in analysis 2, non-parametric kernel, 
bandwidth R = 30 m and 1 m resolution were used. Mutual in- 
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teractions between pairs of species were checked, for example 
species 1 against species 2 (g 12 ) and species 2 against species 1 
(g 2 i) because interactions between species can be asymmetric 
(Wiegand et al. 2007). 

Analysis 4: Species association between life stages 
To investigate spatial associations between life stages of con- 
specific trees, we used the bivariate pair-correlation function 
with null model of random labeling. We considered conspecific 
associations between adults (pattern 1) and individuals of early 
life stages (pattern 2) for each tree species. Parameter g\ 2 -g\\ < 0 
shows a lower intensity of young trees than of adults in the vi¬ 
cinity of adults, while g 2 \-g 22 < 0 indicates that there is additional 
aggregation of young individual independent from conspecific 
adults. 

All point-pattern analyses were performed with the grid-based 
software Programita (Wiegand and Moloney 2004) and the R 
package Spatstat 2.12.2 (Team 2010). Significant departures 
from the null models were evaluated by constructing approxi¬ 
mately 95% confidence envelopes using the fifth lowest and fifth 
highest values of 199 Monte Carlo simulations. 

Results 

Species composition and characteristics 

In plot 1,17 species were found from 12 families, and in plot 2, 
26 species belonging to 14 families. Fourteen species occurred in 
both plots, including four dominant species with more than 50 
individuals (Table 2). 


Table 2: Characteristics of 6 dominant tree species (> 50 individuals/plot) 
in the two 1-ha plots. Babe NP. dbh (diameter at breast height); NN- 
nearest neighbor distance. 


Species 

Density 

(N/ha) 

Basal area 

(m 2 /ha) 

Max dbh 

(cm) 

Mean dbh 

(cm) 

Median 

NN (m) 

Plot 1 






H. kurzii 

246 

10.81 

79.94 

19.27 

2.50 

B. hsienmu 

53 

6.41 

80.57 

35.56 

5.51 

T. ilicifolia 

316 

1.93 

30.57 

7.68 

1.72 

D. sylvatica 

414 

2.95 

25.47 

8.85 

2.25 

S. wightianum 

193 

4.51 

62.74 

14.48 

2.89 

Plot 2 






H. kurzii 

146 

2.23 

46.81 

10.96 

2.69 

B. hsienmu 

101 

5.91 

86.61 

21.64 

2.82 

S. macrophyllus 

397 

4.13 

51.91 

10.09 

1.34 

D. sylvatica 

461 

1.24 

21.65 

7.99 

2.04 

S. wightianum 

219 

1.98 

21.65 

9.60 

2.66 


Two dominant tree species were found in one plot only: T. ilici- 
folia and S. macrophyllus. For analysis 1, all adults of all species 
were used. In all other analyses, we analyzed the six dominant 
species only. In plot 1, H. kurzii had the largest basal area and 
mean dbh. B. hsienmu had the lowest density of the analyzed 
species, but had a large basal area of 6.41 m 2 /ha. T. ilicifolia and 
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D. sylvatica were the most frequent species in plot 1. 

On average, trees in plot 2 were smaller than those in plot 1. 
Densities of S.wightianum and D. sylvatica were similar in both 
plots. For S. wightianum, the considerably larger dbh in plot 1 
resulted in the basal area being twice as large as in plot 2. H. 
kurzzi, D. sylvatica, T. ilicifolia, S. macrophyllus and S. 


wightianum showed a reverse J-shaped distribution in both plots, 
with the number of individuals gradually decreasing with larger 
diameters (Fig. 2). The size distribution of B. hsienmu differed 
between the plots: unimodal in plot 1 (Fig. 2a) and reversed 
J-shaped in plot 2 (Fig. 2b). In both plots, B. hsienmu had no 
individual in size class 60-70 cm. 


(a) dbh distribution in plot 1 [cm] 



o 

I 

% 


S. wightianum 


[0,10] [10,20] [20,30] [30,40] [40,50] [50,60] [60,70] 


300 

250 

200 

150 

100 

50 
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(b) dbh distribution in plot 2 [cm] 
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D. sylvatica 


[0,10] [10,20] [20,30] 


300 

250 
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50 
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T. ilicifolia 


[0,10] [10,20] [20,30] [30,40] 


B. hsienmu 


[0,10] [10,20] [20,30] [30,40] [40,50] [50,60] [60,70] [70,80] 

D. sylvatica ... s - macrophyllus 


[0,10] [10,20] [20,30] 


250 

200 
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100 

50 

0 


[0,10] [10,20] [20,30] 


[0,10] [10,20] [20,30] [30,40] 


Fig. 2: dbh distribution of 6 dominant tree species in the two 1-ha plots, x-axes display size classes in 10 cm intervals, y-axes display number of individ¬ 
ual trees. 
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Environmental heterogeneity trast to this, in plot 2, the /.-function showed a clear departure 

For both plots, we compared the spatial distributions of all adult from CSR for scales r > 5 m (Fig. 4c) and did not approach g = 0. 
trees with dbh >10 cm (Fig. 3) to the null model of CSR to in- In addition, the g -function was significantly > 1 for scales larger 

vestigate if the patterns showed large-scale heterogeneity. Then, than 13-25 m and did not approach g = 1 (Fig. 4d), providing 

the compatibility of Plot 2 to the null model of the HP was evidence of large-scale heterogeneity at plot 2. Moreover, the 

checked with bandwidth R = 30 m. spatial arrangement fit very well with the null hypothesis of HP 

In plot 1, the /.-function showed no large-scale departure from with R = 30 m (Fig. 4e-f) when analyzed by T-function and the 

CSR (Fig. 4a). The g -function fluctuated around g = 1 and fell pair-correlation function. Therefore, for all further spatial-pattern 

completely inside the confidence envelopes for all scales (Fig. analyses, we applied the homogeneous g -function for plot 1 and 

4b). This indicated large-scale homogeneity in the plot. In con- the heterogeneous g-function for plot 2. 


All adults in homogeneous site (plot 1) All adults in heterogeneous site (plot 2) 



0 20 40 60 80 100 0 20 40 60 80 100 

X(m) X (m) 


Fig. 3: Spatial distributions of all adult trees (dbh > 10 cm) in the two 1-ha plots (Size of circles is proportional to dbh of trees). The unit of (x,y) axes is 
meters. 
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Species distributions 

Spatial patterns of dominant species in both sites are shown in 
Fig. 5. All dominant species were found in clustered distribution 
at various scales, except B. hsienmu in the heterogeneous site. 
The aggregation intensity declined clearly with increasing spatial 
scales. All five species were clustered significantly in the homo¬ 
geneous plot (Fig. 5a-e) at various scales under CSR null model. 
No significant pattern was found at scales larger than 30 m. H. 
kurzii, D. sylvatica and S. wightianum were aggregated up to 
large scales (Fig. 5a, c, d). B. hsienmu was random at small 
scales and aggregated at 8-10 m and at larger scales of 20-25 m 


(Fig. 5b). T. ilicifolia showed a strong clustered distribution over 
the entire range of scales up to 27 m (Fig. 5e). 

In the heterogeneous plot, B. hsienmu (inset, Fig. 5i) was the 
only dominant species with random distribution. The four re¬ 
maining dominant species showed aggregated distributions under 
the HP null model with bandwidth R = 30 m (insets, Fig. 5f, j-1). 
In contrast to the homogeneous plot, no significant aggregation 
was found at scales larger than 15 m. This indicates that envi¬ 
ronmental heterogeneity led to different scales in spatial aggre¬ 
gation of a species. 



0 5 10 15 20 25 30 


(k) S. wightianum 




(j) D. sylvatica 



Scale (m) 


Fig. 5: Spatial patterns of dominant tree species shown by the pair correlation function gn with different null models: CSR at the homogeneous site (a-e) 
and HP at the heterogeneous site (insets, f-m). Black lines are observed patterns; grey lines are approximate 95% confidence envelopes. 


Inter-specific associations 

Species associations were analyzed at scales 0-30 m by using the 
bivariate pair-correlation function. In both plots, we perfonned 
40 bivariate point-pattern analyses for all pairs of dominant spe¬ 
cies. Overall, attraction occurred more frequently in the homo¬ 
geneous site, while independence and repulsion predominated in 
the heterogeneous site (appendix 1). 

In the homogeneous site, 10 pairs showed attraction, with 
eight pairs independence and two pairs repulsion (appendix 1). 
Fig. 6 shows examples of species maps and bivariate patterns in 
plot 1. All species pairs had symmetric interactions. T. ilicifolia 
showed negative associations with all other species. D. sylvatica 
and S. wightianum had positive associations with the three re¬ 


maining species. B. hsienmu and H. kurzii were independent in 
species interactions. T. ilicifolia was found only in the lowest 
parts of the plot where it gets higher soil moisture, but is sepa¬ 
rated from other species. In contrast, B. hsienmu and S. 
wightianum are shade intolerant and found in the highest parts of 
the site. 

In the heterogeneous site, four species pairs showed attraction, 
11 pairs independence and five pairs repulsion (appendix 1). Fig. 
7 shows some examples of species maps and bivariate associa¬ 
tion patterns in this plot. Two pairs of species exhibited asym¬ 
metric interactions: B. hsienmu versus S. macrophyllus and ver¬ 
sus S. wightianum, respectively. S. macrophyllus showed only 
negative and independent associations with other species. As in 
the homogeneous plot, S. wightianum had positive associations 
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with B. hsienmu and D. sylvatica. B. hsienmu and H. kurzii 
showed independent associations with the other species, with the 
exception of S. macrophyllus. S. macrophyllus repulsed all spe¬ 
cies except B. hsienmu. 


In the heterogeneous plot, more asymmetric species interactions 
were found. Moreover, the number of significant associations 
was lower and the number of independent species pairs higher 
than in the homogeneous plot. 


(a) T. ilicifolia vs. H. kurzii 

100 7 ojp - 


E, 

>- 



V D .»• 


100 


(c) D. sylvatica vs. H. kurzii 
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50 






ve&il; 

“v> 'tsj 


100 


50 

X (m) 


100 


1.4 

1.2 
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0.6 

0.4 

0.2 

0.0 


100 


(e) S. wightianum vs 8. hsienmu 


E, 

>- 




50 100 

X (m) 

(f) S. wightianum vs. 8. hsienmu 



0.0 


10 15 20 25 

Scale (m) 


30 


Fig. 6: Examples of species distribution maps and bivariate patterns in the homogeneous site. Pattern 1 is displayed as solid circles and pattern 2 as open 
circles. The pair correlation function gn was performed under null model that pattern 1 was fixed and pattern 2 randomized. Black lines are observed 
patterns; grey lines are approximate 95% confidence envelopes. 


(a) S. macrophyllus vs. D. sylvatica 



(c) H. kurzii vs. S. macrophyllus 



50 

X (m) 


100 


0 


100 



50 

X (m) 

(d) H. kurzii vs. S. macrophyllus 


(e) S. wightianum vs. D. sylvatica 





Fig. 7: Examples of species distribution maps and bivariate patterns in the heterogeneous site. Pattern 1 is displayed as solid circles and pattern 2 as open 
circles. The pair correlation function gi 2 was performed under null model that pattern 1 was fixed and pattern 2 randomized under HP process. Black 
lines are observed patterns; grey lines are approximate 95% confidence envelopes. 
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Intra-specific associations between life stages 

Spatial associations between life stages of six dominant tree 
species in both plots were analyzed with the bivariate 
pair-correlation function under null model of random labeling. 
Two test statistics, gn-gn and g 2 i-gn, were used to find con- 
specific associations between adult and younger growth stages 
(appendix 2). 

In the homogeneous plot, B. hsienmu was not analyzed be¬ 
cause of low abundance (less than 30 individuals in each growth 
stage). In H. kurzii and T. ilicifolia, gn-gn < 0 predominated in 
associations between adults and trees at early life stages, indi¬ 
cating that there were fewer young individuals in the vicinity of 
potential parent trees. In D. sylvatica and S. wightianum, gn-gn 
= 0 showed that adults and individuals at early growth stages 
exploited environmental conditions in a similar way. There was 
only one species ( S. wightianum) where juveniles were concen¬ 
trating at higher density around their adults. In H. kurzii, T. ilici¬ 
folia and S. wightianum, g^-gii < 0 dominated, suggesting that 
there was additional clumping of young individual trees inde¬ 
pendent from their potential parent trees. In D. sylvatica, g 2 1 -gn 
= 0 showed no higher density of young individuals in the vicinity 
of theirparents. 

At the environmentally heterogeneous site, only adults and 
sub-adult stages of B. hsienmu were considered, but all life 
stages of the other four species were considered. The test statistic 
gn~g 11 < 0 in all cases of D. sylvatica and S. wightianum indi¬ 
cated a lack of young trees around adult trees, except in D. sylva¬ 
tica where sub-adults concentrated at a higher density around 
conspecific adults (gi 2 -gn > 0). Moreover, g 12 -gn = 0 predomi¬ 
nated in all analyses of H. kurzii, S. macrophyllus and B. 
hsienmu, showing that young and adult trees exploited available 
resources in the same way. Meanwhile, the test statistic g 2 i~g 22 < 
0 in most of the analyzed cases showed that individuals at young 
growth stages were aggregated independently of their conspeci¬ 
fic adults. 

Discussion 

Environmental heterogeneity 

The main goals of our study are to detect effects of environ¬ 
mental heterogeneity on spatial patterns and distributions of tree 
species. For comparative analyses, we use two fully mapped 1-ha 
plots that have four dominant tree species in common, but differ 
in environmental homogeneity. Our results clearly show a dif¬ 
ference in dbh distributions between the individuals of a species 
on the homogeneous and the heterogeneous site. This indicates 
exogenous effects on successful establishment and mortality. 

In the homogeneous site, H.kurzii and S. wightianum present a 
wide range of dbh classes with decreasing numbers of individu¬ 
als in the larger size classes suggesting successful establishment 
and self-thinning in homogeneous site conditions. B. hsienmu 
seems to be well developed in the heterogeneous site with a re¬ 


verse J-shape distribution which contrasts to a lack of individual 
trees in small size classes of the homogeneous plot. That may be 
caused by its ecological preference of direct sunlight as a 
shade-intolerant tree. 

In terns of spatial patterns, the density of all adult trees in the 
homogeneous plot is equal over the entire plot, implying homo¬ 
geneity of environmental conditions. This indicates that the spe¬ 
cies’ spatial distribution depends only on their specific ecology, 
such as seed dispersal or habitat preference. In the heterogeneous 
site, species distribution is affected not only by the species’ 
ecology but also by inhomogeneous site conditions. The inho¬ 
mogeneous pair-correlation function shows the “pure” effect of 
heterogeneity at scales up to the bandwidth R = 30 m. To detect 
the “true” density dependence of trees, it is necessary to separate 
confounding effects. Using adult trees to factor out environ¬ 
mental heterogeneity is a simple and efficient approach, thus it 
has been used in previous studies, for example, Getzin et al. 
(2008) and Zhu et al.( 2009). 

Species distributions 

In our study, all dominant tree species in both plots exhibit ag¬ 
gregated patterns. This is similar to the results of previous stud¬ 
ies, where most tropical tree species were found to have a 
clumped or random distribution (He et al. 1997; Condit et al. 
2000; Bunyavejchewin et al. 2003; Zhang et al. 2013). Spatially 
heterogeneous habitat or limited seed dispersal can lead to 
patchy distribution but may not affect all species (He et al. 1997; 
Harms et al. 2000; Bunyavejchewin et al. 2003; Zhang et al. 
2010 ). 

The univariate analyses of six dominant species reveal sig¬ 
nificant environmental effects on the spatial scale of species 
demographics: H. kurzii, D. sylvatica and S. wightianum are 
aggregated up to larger scales under homogeneous site condi¬ 
tions than under heterogeneous conditions. T. ilicifolia and S. 
macrophyllus are aggregated at large scales indicating their 
habitat preferences and/or well dispersed seeds. Furthermore, T. 
ilicifolia is a humidity-demanding species, favoring relatively 
wet and shady habitat, and is distributed at the lowest area of plot 
1. H. kurzii, S. macrophyllus and D. sylvatica are shade tolerant 
species (FIPI 1996) and they grow under forest canopy. More¬ 
over, D. sylvatica is a dioecious species and can be pollinated by 
wind (Somanathan and Borges 2000) which can explain its 
strong aggregation. B. hsienmu is mainly dispersed by gravity 
and has a poor regeneration (Wang et al. 1986) leading to a low 
number of young trees (plot 1) in a random pattern (plot 2). 

Light conditions are sufficient for S. wightianum in canopy 
gaps since it is a light-demanding species (FIPI 1996; 
Gunatilleke et al. 2006). Direct intra-specific interactions occur 
within scales of 30 m and disappear at larger scales (shown by 
Plot 1). This finding is consistent with previous studies using 
individual-based analyses, where local neighborhood effects 
occurred within 20 or 30 m (Hubbell et al. 2001; Uriarte et al. 
2004b; Wiegand et al. 2007). 
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Inter-specific associations 

Attraction between species predominates in the homogeneous 
site while repulsion is more frequent in the heterogeneous site. 
Positive association suggests that species have similar habitat 
preference or occupy different forest layers. The latter leads to 
different regeneration requirements (Yamada et al. 2006). The 
attractions between B. hsienmu - S. wightianum (plot 1, both 
light demanding but B. hsienmu is a slow-growing species), D. 
sylvatica - S. wightianum (in both plots, D. sylvatica is shade 
tolerant and a sub-canopy tree) and B. hsienmu - D. sylvatica can 
be explained by their ecological differences. Negative interac¬ 
tions emphasize the effect of habitat segregation onto species 
association in forest stands. Among the three main limiting plant 
resources, soil nutrients, water availability and light (Debski et al. 
2002 ), we assume that soil humidity and light are the driving 
factors in the study sites. 

Segregated distribution or topographical habitat differentiation 
can be reasons for negative interactions (North et al. 2004; 
Pizano et al. 2011). This could explain why T. ilicifolia (plot 1) 
and S. macrophyllus (plot 2), whose distributions are separated, 
show negative associations with most other species. 
Bunyavejchewin et al.(2003) found that negative associations 
provide indirect evidence of habitat segregation. H. kurzii shows 
differences between the plots regarding its associations with 
other species: it is attracted by other species in the homogeneous 
site while it is independent from other species in the heteroge¬ 
neous site. Habitat segregation is based on the niche differentia¬ 
tion of resources; therefore a segregated habitat could be suited 
for different species and leads to less inter-specific associations 
(Harms et al. 2001). The fact that independence dominates in 
species associations at the heterogeneous site might be due to the 
effects of biotic and abiotic factors. Assertion of species inde¬ 
pendence may be a good approximation for species-rich commu¬ 
nities, but still remains unclear in communities with low species 
richness (Wiegand et al. 2012). 

Intra-specific associations between life stages 

Bivariate analyses between conspecific life stages provide a 
closer view on intra-specific associations. In some species, 
gn-gn < 0 implies a lack of young individuals around conspeci¬ 
fic adults, such as H. kurzii, T. ilicifolia (plot 1) and D. sylvatica, 
S. wightianum (plot 2). Moreover, gi 2 -gn — 0 in the remaining 
species indicates that young trees surround adult trees in the 
same way as adults surround adults. Therefore, these results 
imply low densities of “recruits” in the vicinity of their conspeci¬ 
fic potential parents. This lack of recruits near conspecific adults 
might be consistent with Janzen-Connell effects at small scales 
(Wright 2002) and establishment limitation due to abiotic or 
biotic constraints (Luo et al. 2011). The Janzen-Connell hy¬ 
pothesis describes host specific pests and or pathogens reduce 
recruitment near conspecific adults (Janzen 1970; Connell 1971). 
On the other hand, seeds that get far from parent trees have more 
chance to survive and arrive in favorable sites (Hood et al. 2004). 
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Additional support for the Janzen-Connell hypothesis is provided 
by the fact that g 2 rg 22 < 0 for all life stage associations, except 
those in D. sylvatica implying that young individuals are 
clumped independently and not closely to their conspecific adults 
and highlighting recruitment limitation near conspecific adults. 

In conclusion, our major objectives were to explore the spatial 
patterns and associations of dominant tree species and how they 
are influenced by heterogeneous environmental conditions. 
Through spatial point-pattern analyses, the results clearly show 
the impacts of environmental heterogeneity on horizontal struc¬ 
ture, and intra- and inter-specific interactions. Meanwhile, there 
are still many unanswered questions that may be solved if more 
infonnation is integrated on the ecological characteristics of trees. 
Moreover, spatial pattern analysis at larger scales may exhibit 
more effects of environmental heterogeneity as suggested by 
Boyden et al.(2005). 
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Appendix 1: Inter-specific association: (+) attraction, (0) independence and (-) repulsion 
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Appendix 2: Intra-specific association between life stages: (+) attraction, (0) independence and (-) repulsion 
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